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THE EFFECTS OF CHARGE TRANSFER INEFFICIENCY (CTI) ON GALAXY SHAPE MEASUREMENTS 

Jason Rhodes 1,2 , Alexie Leauthaud 3,4 , Chris Stoughton 5 , Richard Massey 6 , Kyle Dawson 7 , William Kolbe 3 , 

Natalie Roe 3 

ABSTRACT 

We examine the effects of charge transfer inefficiency (CTI) during CCD readout on the demanding 
galaxy shape measurements required by studies of weak gravitational lensing. We simulate a CCD 
readout with CTI such as that caused by charged particle radiation damage in space-based detectors. 
We verify our simulations on real data from fully-depleted p-channel CCDs that have been deliberately 
irradiated in a laboratory. We show that only charge traps with time constants of the same order as 
the time between row transfers during readout affect galaxy shape measurements. We simulate deep 
astronomical images and the process of CCD readout, characterizing the effects of CTI on various 
galaxy populations. Our code and methods are general and can be applied to any CCDs, once the 
density and characteristic release times of their charge trap species are known. We baseline our study 
around p-channel CCDs that have been shown to have charge transfer efficiency up to an order of 
magnitude better than several models of n-channel CCDs designed for space applications. We predict 
that for galaxies furthest from the readout registers, bias in the measurement of galaxy shapes, Ae, 
will increase at a rate of (2.65 ± 0.02) x 10~ 4 yr~ x at L2 for accumulated radiation exposure averaged 
over the solar cycle. If uncorrected, this will consume the entire shape measurement error budget 
of a dark energy mission surveying the entire extragalactic sky within about 4 years of accumulated 
radiation damage. However, software mitigation techniques demonstrated elsewhere can reduce this 
by a factor of ~ 10, bringing the effect well below mission requirements. This conclusion is valid only 
for the p-channel CCDs we have modeled; CCDs with higher CTI will fare worse and may not meet 
the requirements of future dark energy missions. We also discuss additional ways in which hardware 
could be designed to further minimize the impact of CTI. 

Subject headings: cosmology: observations - gravitational lensing - large-scale structure of Universe 



1. INTRODUCTION 

The past decade has seen enormous changes in the field 
of cosmology. A concordance cosmology in which the ex- 
pansion of the universe is accelerating has been accepted 
(Spergel et al. 2007). This accelerated expansion was 
first demonstrated by observations of SN la (Riess et al. 
1998; Perlmutter et al. 1998) and has been confirmed 
with other probes in recent years. This startling aspect 
of our Universe has prompted a wide variety of possi- 
ble explanations (see, e.g., Caldwell 2004) and consid- 
erable effort has gone into developing concepts for dedi- 
cated space missions to probe the mysterious dark energy 
thought to be causing this accelerated expansion. These 
missions, which include the NASA/DOE Joint Dark En- 
ergy Mission^ and ESA's Euclid mission^, plan to use a 
variety of probes in order to constrain the properties of 
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dark energy. It has become widely accepted that weak 
gravitational lensing, the small distortion of the observed 
shapes of background galaxies by foreground dark mat- 
ter, is one of the most powerful probes of dark energy, 
provided that systematic effects can be controlled (Al- 
brecht et al. 2006). Space missions are attractive in large 
part due the greater ability to control many systematic 
effects (Rhodes et al. 2004a). 

The field of weak lensing by large-scale structure, or 
cosmic shear, developed in parallel to the study of dark 
energy as the dominant component of the Universe. 
From the first detections of cosmic shear a decade ago 
(Wittman et al. 2000; Bacon, Refregier & Ellis, 2000; 
Kaiser, Wilson, & Luppino, 2000; Van Waerbeke et al. 
2000), surveys have grown in size and thus information 
content (see Hoekstra & Jain 2008 for a recent review). 
The culmination of this effort will be in the execution 
of the above-mentioned dedicated dark energy missions, 
which plan to survey up to 20000 square degrees, the 
entire extragalactic sky. It has become clear that con- 
trol of systematic effects, both observational and astro- 
physical, is of paramount importance in making use of 
weak lensing as a probe of dark energy. The subtle shape 
changes induced by weak lensing require exquisite control 
of observational systematic effects, especially knowledge 
of the telescope's point spread function (PSF). From the 
ground, thermal and gravity load-induced fluctuations in 
the telescope can change the PSF, and the atmospheric 
seeing both broadens the PSF and makes the PSF un- 
stable on the timescale of astronomical exposures. These 
effects can be largely or completely mitigated by making 
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observations in a thermally stable environment above the 
atmosphere. 

The control of systematics is a large factor in the drive 
towards a dedicated space mission. However, the harsh 
radiation environment of space has the potential to in- 
troduce an observational systematic effect due to charge 
traps created in CCD detectors by impacts from charged 
particles. These defects trap charge (either electrons in 
so-called n-channel CCDs or holes in p-channel CCDs) as 
charge is clocked across pixels toward the readout regis- 
ters. When the charge is subsequently released from the 
trap, it shows up in a neighboring pixel, thus creating 
a trail along the readout direction. These trails obvi- 
ously change the observed shapes of the galaxies in the 
images. These shape changes are coherent across the im- 
age, thus mimicking a weak lensing signal caused by dark 
matter. The degradation of charge transfer efficiency 
(CTE, this quantity is one minus the CTI, or charge 
transfer inefficiency) due to these radiation-damage in- 
duced traps has been observed in all Hubble Space Tele- 
scope (HST) imaging cameras: the Wide Field Planetary 
Camera 2 (WFPC2; Dolphin 2009), the Space Telescope 
Imaging Spectrograph (STIS; Goudfrooij et al. 2006), 
and the Advanced Camera For Surveys (ACS; Sirianni 
et al. 2005) and has hampered the sensitive shape mea- 
surements needed for weak lensing with those cameras 
(Rhodes et al. 2004b; Schrabback et al. 2007; Rhodes 
et al. 2007). CTE degradation is a particularly difficult 
effect to correct for because its non-linear nature means 
that high signal-to- noise (S/N) stars which are typically 
used in weak lensing for PSF modeling will be affected 
less than the low S/N galaxies whose shapes are being 
measured. Thus, typical PSF deconvolution techniques 
are complicated by the effects of CTI. Thus, it is clear 
that future space weak lensing missions will need to min- 
imize CTI due to radiation damage and have CCDs that 
are sufficiently well understood to allow for mitigation of 
the CTE degradation that does occur. 

In this paper we carry out a quantitative analysis of the 
effects of CTI in CCDs on the analysis of galaxy shapes 
for weak lensing, and explore techniques to mitigate the 
shape distortions due to trailing charge. We have de- 
veloped a detailed model for the effect of charge traps 
based on data from irradiated CCDs, and applied it to 
simulated galaxies. We quantify the effects on measured 
galaxy shapes in simulated data as a function of galaxy 
size, CTI, and S/N. We base our analysis on data from p- 
channel CCDs fabricated at Lawrence Berkeley National 
Laboratory (LBNL) and irradiated with protons at the 
LBNL 88" cyclotron (Dawson et al 2008; hereafter D08). 
However, our methods are general and can be applied to 
any CCDs if the density and time constants of the charge 
traps are known. 

This paper is organized as follows. $2] gives a brief 
overview of how charge is transferred between pixels dur- 
ing CCD readout and how different types of CCDs have 
performed after radiation damage. £j3] describes the code 
we have developed to mimic the effects of CTI on CCD 
readout. In S|4]we validate that code by showing that it 
can reproduce the effects of CTE degradation as mea- 
sured on real, irradiated LBNL CCDs. We apply our 
code to simulated astronomical images and detail the ef- 
fects of CTI on the shapes of galaxies in fj5] We examine 
how this will effect the future space missions and give rec- 



ommendations for mitigating the effects in [JSJ Finally, 
we offer concluding remarks in fj7] 

2. BACKGROUND 
2.1. CCD readout and charge transfer 

During exposure, a photon incident on a CCD gener- 
ates an electron-hole pair. In n-channel CCDs, the elec- 
trons drift into the potential well of the nearest pixel, 
which is created by an electrostatic potential gradient 
within the substrate. In p-channel CCDs, holes are col- 
lected instead, but we shall not distinguish between the 
two mechanisms hereafter. Charge (electrons or holes) 
accumulates in a well-defined volume, outside which the 
density falls rapidly to zero (Hardy, Murowinski & Deen 
1998, Seabroke et al. 2008). We are most concerned with 
the cross-sectional area of this cloud, which expands as a 
monotonic function of the amount of charge n e . As illus- 
trated in Figure [TJ we parameterize this as an effective 
height h within the pixel; the variable need not really be 
the height, but that is a useful one-dimensional function 
for the purposes of explanation. 

Some CCDs contain a supplementary buried channel 
or "notch" constructed at the bottom of the potential 
well. This is a small region of slightly lower potential 
and, like the channel in an artificial river bed designed 
to improve the flow of a small amount of water, it will 
concentrate the first few electrons (or holes). The notch 
reduces the number of traps the charge packet is exposed 
to, by confining the charge to a smaller region. We model 
the notch by setting the height h to zero below some 
notch depth d, and 

above the notch, where w is the full well depth (the total 
amount of charge that will fit in a pixel before it over- 
flows) and a, which depends on the construction of the 
potential but which is typically ^0.5 (see for instance, 
Chiaberge et al. 2009 and Mutchler & Sirianni, 2005). 
Setting h to zero below the notch is an approximation 
because it is possible that there are some traps in the 
notch. At any rate, the CCDs we model in this paper 
have no notch in the imaging region and thus this ap- 
proximation does not affect the results presented here. 

At the end of the exposure, the charge must be moved 
to the amplifier electronics at the edge of the CCD to be 
read out and digitized. Each row of charge is first shuf- 
fled one pixel in the parallel direction, towards the serial 
readout register. This is typically accomplished using a 
3-phase clock, as is the case for the CCDs described here 
(Janesick 2001), but can be approximated as a single 
operation (for a discussion of the consequences of this 
approximation, see Massey et al. 2009). Charge from 
the bottom-most row is transferred into the serial read- 
out register and then shuffled using the same technique 
but in a perpendicular direction. The charge from each 
pixel is shifted onto a capacitor connected to an ampli- 
fier and then 'counted' by being sensed as a voltage and 
digitized. Read noise is the shot noise on this voltage 
and can be reduced by lengthening the sampling time, 
and thus slowing the readout rate. This process is then 
repeated once for each additional row of pixels: shuffling 
one row in the parallel direction and then through the 
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serial readout register. 

As the charge is transferred from pixel to pixel, charge 
traps due to defects or impurities in the Si lattice tem- 
porarily capture passing charge, and release it after some 
delay. The typical capture time is effectively instanta- 
neous (and if it is not, the lower capture rate can be 
equivalently modeled as a lower density of charge traps) . 
The probability of release is governed by an exponential 
decay, with a characteristic time constant that depends 
upon the properties of the lattice and impurities, as well 
as the operating temperature of the detector (Shocklcy 
& Read 1952, Hall 1952). Several different species of 
charge traps may be present in any given device, with 
different characteristic release times r, and densities p. 
If the captured charge is released after its charge cloud 
has been shuffled along, the released charge appears as a 
faint trail behind the main charge packet. 

Trapping and trailing can occur in both the parallel 
and serial directions. Each column of pixels is indepen- 
dent and has a unique set of charge traps, while all of 
the rows share the same traps along the serial register. 
Note that even though the pixels and transfer mecha- 
nisms are physically similar in the two directions, the 
time per pixel transfer is typically ~ 10 3 shorter in the 
serial direction (D08). We demonstrate in 52] that only 
charge traps with release times roughly similar to the 
clocking (time between transfers) affect galaxy shapes. 
Thus, separate species of charge traps can be important 
for parallel and serial transfers. 

2.2. CTE Effects in Irradiated CCDs 

Thick, fully depleted p-channel CCDs have several ad- 
vantages over conventional thin, n-channel CCDs, includ- 
ing enhanced quantum efficiency at near-infrared wave- 
lengths, reduced fringing at near-infrared wavelengths, 
and significantly less degradation of CTE with a given 
accumulated radiation exposure. This last effect is par- 
ticulary important for the measurement of subtle weak 
lensing-induced galaxy shape distortions. It arises due to 
the fact that the divacancy traps that are primarily re- 
sponsible for CTI in p-channel devices are more difficult 
to form than the phosphorous vacancy traps that occur 
due to radiation damage in n-channel CCDs (Bebek et al. 
2002; Janesick & Elliott 1992; D08; Spratt et al. 2005). 
We focus our analysis on p-channel CCDs precisely be- 
cause it is more difficult to form traps in these types of 
CCDs than in existing n-channel CCDs. 

Marshall et al. (2004) compared the CTE responses of 
irradiated p-channel LBNL CCDs with n-channel CCDs. 
These n-channel devices are designed for space applica- 
tions and are the ones used for the recently installed 
Wide Field Camera 3 on the HST. They found that a 
notch implant in the channels improved the CTE perfor- 
mance by a factor of 2 for both p-channel and n-channel 
devices. More importantly, they found that the CTE 
performance of p-channel devices is about an order of 
magnitude better than that of n-channel devices after 
irradiation. A re-analysis of the Marshall et al. data 
(Lumb 2009) indicates that the p-channel devices may 
only have about a factor of 3-8 better CTE (depending 
on the signal level, with the p-channel advantage being 
greater at low signal levels such as those expected in the 
images of faint galaxies). Likewise, a comparison by Gow 
et al.(2009) found a factor of 7 improvement in tolerance 



for parallel CTI (and 15 for serial CTI) in otherwise sim- 
ilar p and n-channel devices. Thus, thick, fully depleted 
p-channel CCDs are particularly attractive for a weak 
lensing space mission and we use them as the baseline 
for this study. We do, however, note that p-channel de- 
vices do not have the rich heritage that n-channel devices 
do, particularly in space applications. 

LBNL has developed radiation-hardened CCDs with 
the specific application of dark energy missions in mind 
(Holland et al. 2006). These CCDS are composed of 
3512 x 3512 10.5/im pixels with 4 readout registers and 
were baselined for the SuperNova Acceleration Probe 
(SNAP), a JDEM concept (Bebek 2007). The SNAP mis- 
sion, like other candidate dark energy missions, would be 
at the L2 Earth-Sun Lagrange point, and we use the ra- 
diation flux there as the baseline for the flux experienced 
by a dark energy mission. A dark energy observatory 
will experience significant radiation exposure at L2, pri- 
marily from solar protons. Exposure to energetic protons 
leads to degraded CCD performance due to bulk dam- 
age from non-ionizing energy loss (NIEL) and charging 
of oxide layers from ionizing radiation. Bulk damage in 
the Si lattice is the dominant effect that manifests itself 
through increased CTI, increased dark current, and iso- 
lated hot pixels. Of these effects, an increase in CTI is 
the most likely to introduce systematic errors to a weak 
lensing survey. 

3. SOFTWARE ALGORITHM TO MIMIC CTI 

Our model for CTI is that inefficient charge transfers 
are caused by discrete charge traps embedded in pix- 
els. These charge traps can capture and release single 
electrons (or holes). Each pixel of the detector contains 
a number of charge traps, that have a mean density p. 
Each trap is characterized by htrap, its vertical location 
in the pixel and r, its characteristic release time con- 
stant. Different species of charge traps have different 
decay time constants and thus different values of r (and 
p). Whenever a trap in a pixel containing n e electrons 
(or holes) is within the electron cloud, i.e. htrap < h(n e ), 
we assume that it immediately absorbs one electron from 
the free charge. When a full charge trap is above the free 
charge height (that is, not within the electron cloud), it 
is given an opportunity to decay. 

We have developed a code that mimics the readout of 
a CCD with imperfect CTI. For each column of pixels, 
we use the following procedure to "read out" the image 
and determine the observed charge in each pixel. This 
procedure is also illustrated graphically in Figure [2] 

1. Populate each pixel with free charge, as accumu- 
lated during an exposure, and calculate h(n e ). 

2. Define the locations of charge traps throughout the 
pixel array. This is done for multiple species of 
charge traps, each with a different p and t. 

3. For charge traps with htrap < h(n e ) fill the trap 
with one electron, and subtract one electron from 
the free charge in its pixel. 

4. Read out the charge in pixel row n = 1. 

5. For n = 1 to n = n ma x set the free charge in pixel 
n to the free charge in pixel n + 1 . 
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Fig. 1. — Cartoon illustrating the CCD readout process. 



1. Populate each pixel 
with (free) charge 



2. Define locations 
of charge traps 



3. Capture any free 
charge below the 
height of a trap 



6. Probabilistically 
release trapped charge 
above height of traps 



4. Read I 
out serial 
register I 



5. Move all free charge 
to adjacent row 



Y 

Finish 

Fig. 2. — Flowchart representing the readout process as imple- 
mented in our code. 

6. For full traps with h trap > h(n e ), calculate the 
probability that the trap will decay (that is, release 
the charge), based on an exponential decay with 
time constant r. Generate a random number in 
the range [0,1] and, if the probability is less than 
the random number chosen, empty the trap and 
increase the free charge by 1 electron. 

7. Repeat the previous four steps to calculate the 
measured charge in pixel rows n = 2, n — 3, .... 
Note that n max is decreased by one for each itera- 
tion. 

The image is then rotated by ninety degrees and ex- 
actly the same process is repeated to simulate serial 



transfers. The total number of operations to read out 
one row or column of n P i X pixels is (n P i X )(n P i X + l)/2. 
Thus, reading out out an entire n P i X x n P i X array scales 
as rip ix and is computationally intensive. During a read- 
out, a charge can be trapped multiple times in different 
pixels. 

This process inevitably adds noise to an image because 
we do not know the true locations of individual traps 
within the detectoiFl. but simply model the traps as a 
uniform density. We therefore ought to perform many 
iterations and average the results. As a more compu- 
tationally efficient solution to this problem, we instead 
introduce "fractional traps" . That is, we place the same 
number of traps in each pixel and let each trap capture a 
"fractional electron" (with the fraction being p, the trap 
density per pixel, which can be less than 1). We fur- 
ther divide the fractional trap in each pixel into ni eve i s 
fractional traps, each located at a vertical position in 
multiples of l/ni eve i s . These traps release charge ex- 
actly as described above. We have found that setting 
ni eve i s — 10000 allows us to reproduce the averaged re- 
sults of many iterations with full traps placed randomly 
within the pixels (mimicking a real, physical CCD). How- 
ever, the fractional trap method saves considerable com- 
putational overhead when simulating the effects of CTI 
on many thousands of images as described in Sj5j Our 
code allows each trap to have a different t. However, 
as discussed below, we find that in each charge trans- 
fer direction, a small number of t values (trap species) 
describe the behavior of physical detectors. 

4. VALIDATION OF THE CHARGE TRANSFER CODE 

The charge transfer code has been tested and validated 
by creating simulated images designed to mimic the ir- 
radiated data used in D08. Using the same software as 
D08, we show that the code is able to reproduce the ob- 
served radiation damage effects as measured in the D08 
data. 

4.1. Irradiation in the Lab to simulate Space Radiation 

10 This knowledge may be possible in future analyses using 
"pocket pumping", as described in [0 
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The damaging particle radiation flux incident on a 
spacecraft depends on exactly when the mission occurs 
within the ~ 11 year solar cycle (Barth et al. 2000). 
There is an order of magnitude difference in the flux of 
solar protons during the heaviest and lightest parts of 
the solar cycle (see, e.g. Figure 1 of Barth et al. 2000). 
In D08, the solar proton flux was modeled using the Eu- 
ropean Space Agenc y 's Space Environment Information 
System (SPENVIS£3- In SPENVIS, a simplified solar cy- 
cle consisting of 7 years near the maximum flux level and 
4 years at zero flux is used. The total displacement dam- 
age (energy deposited in the silicon) is predicted from a 
uniform 4-7T steradian spatial distribution of solar protons 
with an energy distribution derived from the Xapsos et 
al. (1999) model for solar proton emission. SPENVIS 
employs a statistical model based on data from previ- 
ous solar cycles to predict the dose at 95% confidence 
level (CL); that is, the prediction will underestimate the 
dose only 5% of the time. D08 used the SPENVIS model 
to calculate the 95% CL solar proton flux incident on 
the CCDs, after passing through the shielding provided 
by the SNAP spacecraft and telescope, yielding an inte- 
grated NIEL exposure of 2.54 x 10 6 MeV/g (Si) for one 
year at solar maximum. 

D08 then characterized the CTE performance of thick, 
fully depleted LBNL p-channel CCDs by irradiating sev- 
eral CCDs at the LBNL 88- Inch Cyclotron with 12.5 and 
55 MeV protons. Although a variety of irradiation levels 
were used, we only consider the 12.5 MeV data with an 
irradiation level of 2 x 10 10 protons/cm 2 . The NIEL fac- 
tor for 12.5 MeV protons is 8.9 x 10 _3 MeVg _1 cm~ 2 per 
proton (Jun et al. 2003). Thus, the 2 x 10 10 protons/cm 2 
flux of 12.5 MeV protons used to irradiate the CCDs cor- 
responds to 1.78 x 10 10 MeV/g (Si,) a total dose equiv- 
alent to ten solar cycles at 95% CL, using the SPENVIS 
approximation detailed here, or 110 years at L2. 

Although this accumulated radiation exposure is sig- 
nificantly higher than any proposed dark energy mis- 
sion would encounter, the exaggerated radiation expo- 
sure makes it easier to characterize the detailed effects of 
CTE degradation. Since the number of traps (and thus 
the degradation of CTE) is linear with radiation expo- 
sure and the NIEL dose, the D08 accumulated radiation 
exposure can be used to estimate CCD performance over 
the course of a dark energy mission lifetime. 

4.2. Analysis of CTI Due to Radiation Damage 

The CTE of irradiated CCDs was measured using a 
55 Fe x-ray source that emits K- alpha photons with en- 
ergies of 5.9KeV. At the operating temperature of 133K, 
a single K-alpha x-ray will generate 1580 electron/hole 
pairs, which, depending on the location of the x-ray rel- 
ative to the pixel potential wells, may be localized in a 
single pixel or shared among two or more pixels. 

In D08, CTE was characterized using single pixel 
events from the K-alpha peak, and the results showed 
that the irradiated LBNL CCDs are three times more af- 
fected by charge trailing in the parallel readout direction 
that in the serial readout direction. In this article, we 
disregard the serial CTE and only consider the trailing 
in the parallel direction since this will most affect galaxy 
shape measurements. 

11 See http://www.spenvis.oma.be 



The effects of irradiation on CTE was studied in two 
ways in D08. In the first method, called the stacking 
method, CTE is characterized by the average charge col- 
lected for single pixel x-ray events as a function of the 
number of pixel transfers. Those x-ray events that ex- 
perience more transfers lose a larger amount of charge 
due to CTI, as shown in Figure [3] The serial and par- 
allel CTE components are determined independently by 
fitting the fractional loss of each transfer to the data. 
The same single pixel x-ray events are used for a mea- 
surement of CTE using the trailing method, in which 
the charge is counted in each trailing pixel as a frac- 
tion of the charge in the primary charge packet. The 
fractional trailing charge in each event is divided by the 
total number of transfers, and the results are averaged 
over all x-ray events. In other words, the averaged trails 
represent the fraction of charge left behind the primary 
charge packet for a single transfer. The effect of this 
trailing with a best fit to the data after irradiation and 
1650 transfers is shown in Figure HJ The total fractional 
charge integrated over these trails represents the CTI. To 
summarize and compare the two methods, the stacking 
method offers a direct measurement of CTE, measuring 
charge that is successfully transferred relative to the ex- 
pected x-ray charge deposition, while the trails method 
offers a direct measurement of CTI by measuring the 
trailing charge relative to the charge in the leading pixel. 
The total fractional charge integrated over these trails 
represents the CTI. 

In practice, the analysis of the trails following x-ray 
events is limited by the ability to measure the faint trails 
after a large number of transfers in the presence of non- 
zero read noise. Because of the low S/N at large dis- 
tances from the primary x-ray event, D08 fit only the 
first 45 pixels with a two term exponential. The charge in 
this two-term exponential represents approximately 2/3 
of the total charge lost due to CTE effects as identified in 
the stacking method. In a re-analysis of the D08 data, in 
which stacking plots were made with varying x-ray flux, 
we find strong evidence that the remaining charge must 
be attributed to one or more populations of traps with 
a much longer time constant. One such candidate is the 
C-0 trap which has been independently identified in the 
LBNL CCDs in a previous analysis (Bebek et al. 2002) 
with a de-trapping time constant of many seconds com- 
pared to the typical time between pixel row transfers of 
25 ms at the 70, 000 pixel/sec readout speed employed in 
D08. 

We therefore choose to model the CTI using three dis- 
tinct trap populations instead of the two used by D08. 
We assign the third trap a time constant corresponding 
to 200 pixels in our analysis. The results of our fit are 
shown in Fig. [4] and the best- fit parameters are found to 
be: pi = 0.35, p 2 = 0.49, p 3 = 0.7, n = 10.0, r 2 = 0.486, 
r 3 = 200.0, where r is in units of 'pixels'. These best fit 
parameters are tabulated in Table [T] There are, on av- 
erage, about 1.5 traps per pixel when summed over the 
three species; a mission that was at L2 for only about 
5 years would thus have only about one twentieth that 
number of traps. We created a set of simulated images 
with the same characteristics of the D08 images to test 
our models. Overall, the agreement is good, as shown in 
Figure HJ 

As we show in [J5l charge traps with time constants that 
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TABLE 1 
Parameters of Charge Traps 



Trap Species 


p (traps per pixel) 


t (pixels) 


r (ms) 


probable delect type 


1 


0.35 


10.0 


250 


Carbon-interstial (Ci) 


2 


0.49 


0.486 


12 


divacancy (VV) 


3 


0.7 


200.0 


5000 


carbon-oxygen (CO) 



Xray energy (in units of electrons) 




parallel CTI (stacking) = 6. 12140e-05 +/- 1.64319e-07 
parallel CTE (stacking) = 0.999939+/- 1.64319e-07 



200 400 600 800 1000 1200 

Number of transfers 

Fig. 3. — Decrease in energy of single pixel events as a function 
of the number of transfers in the simulated images. The magenta 
solid line shows the mean energy of the simulated K-alpha impacts 
(1580 e~ ). Because of charge trailing, the measured energy of these 
events after readout decreases with the number of transfers. The 
best fit line to the data is shown by the blue solid line and yields 
a CTE of 0.999939. 

are long relative to the time between parallel transfers in 
readout have negligible effect on galaxy shapes (i.e. they 
do not leave trails). The trapped charge is, however, 
removed from the object, and will thus affect photometry. 
Because our primary motivation is to study the impact 
of charge trailing on galaxy shapes, we chose to mimic 
the measured trails in D08 instead of matching the CTE 
as measured by the stacking method. We note that the 
opposite approach would be appropriate if our aim was 
to measure changes in photometry instead of shapes. 

5. APPLICATION OF THE CTE CODE TO SIMULATED 
GALAXY IMAGES 

We create simulated galaxy images with de Vau- 
couleurs and exponential profiles and use the code de- 
scribed in ij3]to introduce the effects of CTI on the galaxy 
images. We create galaxies as they would appear in the 
I-band of the proposed SNAP mission, which has a 2 
meter mirror, 0.1" pixels and 400 second exposures. The 
background level is chosen to be the average background 
for extragalactic observations taken from L2. The mea- 
sured background is slightly lower than the input back- 
ground because the CTI causes the flux to be dragged 
out into the overscan region of the CCD during readout. 
We create single galaxies in each image to avoid having 
the traps in a pixel be filled by charge from an object 
that has already passed through that pixel during read- 
out. All objects are placed 1650 pixels from the readout 
register (close to the maximum of 1712 pixels from the 
readout register). This is done because objects farthest 
from the readout registers will encounter the most traps 
and suffer the worst CTI; we are trying to estimate the 
worst case scenario for galaxy shape measurement. 

We assume that the trap densities p increase linearly 
over time, as traps accumulate due to radiation damage, 



and that the time release constants r do not change be- 
cause they are properties of the detector material itself. 
The assumption of linearity in trap density is not entirely 
accurate because the proton flux is dominated by solar 
radiation, which varies over the solar cycle. A specific 
analysis of any future dark energy mission will need to 
take into account the portion of the solar cycle in which 
the mission occurs. We make a further simplification 
by assuming that at the start of a mission a CCD will 
have no charge traps, so only traps accumulated during 
the time spent at L2 affect readout. This is of course 
not true, because real CCDs always have some imper- 
fections even immediately after their production. In the 
case of the LBNL CCDs we are simulating, however, this 
turns out to be a good approximation because the pre- 
irradiation CTE is so high and the number of traps so 
low (see Table IV of D08). However, our results repre- 
sent a best-case scenario for the number of traps (and 
thus CTI) as a function of time; the real CTI will be 
slightly worse. 

For each galaxy, we measure the shape both before and 
after the image is degraded with imperfect CTE. The 
shape is parameterized in the typical weak lensing fash- 
ion by a two component ellipticity where ei = y 2 - t \ vv 
corresponds to elongation along the x axis (for positive 



21,. 



cor- 



ei) or the y axis (for negative ei), and e-i 

responds to elongation at ±45 degrees. Here, the shapes 
are described in terms of the second order moments of 
the pixel intensity / such that Jy = ^f^ 3 where Xi is 

the distance in pixels from the object centroid and w is 
a Gaussian weight function. The ellipticities were mea- 
sured using the method of Rhodes, Rcfrcgicr, & Groth 
(2000; hereafter RRG). This method has been well-tested 
on real and simulated space-based data (see Leauthaud et 
al. 2007). Since we are only interested in perturbations to 
the galaxy shapes, we do not go through the somewhat 
complicated steps of point spread function deconvolu- 
tion, which can introduce biases in shape measurements 
(Heymans et al. 2006; Massey et al. 2007). Instead, we 
only concern ourselves with Ae, which is relatively inde- 
pendent of the particularly shear measurement method. 
We explore the effects of S/N and galaxy size on Ae 
but for the bulk of our analysis, we consider small, faint 
galaxies; any weak lensing survey will be dominated by 
galaxies that are faint and small relative to the PSF size. 

For the purposes of this paper, we only introduce par- 
allel CTI into the simulated images and we set the serial 
CTE equal to 100%. We do this for two reasons. First, 
the parallel CTI smears objects in the vertical direction 
(negative ei), but serial CTI smears them horizontally 
(positive ei, the serial readout direction). Thus, just us- 
ing the ellipticity e as an indicator of CTI-induced galaxy 
shape changes means that the effects of serial and paral- 
lel CTI partially cancel. The combined effects change the 
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Fig. 4. — This plot shows the trailing pixel charge after 1650 transfers relative to the charge in the initial charge packet (the first pixel 
of an x-ray event). Left panel: Charge trailing measured behind single pixel events in the parallel read-out direction from DOS (combined 
data set from 14-17 August 2006) (magenta diamonds). The trails in this plot have been calculated up to 100 pixels after the main x-ray 
packet as opposed to 45 pixels in Figure 5 of D08. The dashed black line shows the two trap model used in D08 that underestimates the 
charge trapped in the far reaches of the tail. In order to best fit the data over the range to 100 pixels, we use a 3 population model and 
we determine the trap densities and time constants that best match the data. The predicted trail from our model is shown by the blue 
solid line. Right panel: We create a series of images designed to mimic the D08 data in terms of the background noise and the K-alpha 
impacts (signal). We read out the simulated data with imperfect CTE and measure the charge trailing using the same software as DOS. 
The trails measured in the simulated data are shown by the dark blue asterisk data points and are a good match to the real data (magenta 
diamonds). The error bars for the simulated data would be about the same size as for the real data; for clarity the simulated error bars 
are not plotted. 



size of the PSF, and thus must be corrected for in real 
images, but would provide an unfair test for these pur- 
poses. The second reason we concentrate just on parallel 
CTI is that the parallel CTE degradation is three times 
worse for a given radiation exposure (D08) and different 
charge trap species affect the parallel and serial CTE be- 
cause of the different clocking times in the parallel and 
serial directions. Thus, we seek only to demonstrate that 
we understand the more influential effects of the parallel 
CTI on shape measurement in this paper. 

5.1. Effect of charge trap release time 

Figure [5] demonstrates the effect of charge trap release 
time t on the measurement of photometry as measured 
by S/N (calculated via SExtractor; Bertin & Arnouts 
1996, S/N =flux_auto/fluxerr_auto), astrometry (the y 
centroid of an object), and shapes (rms size d rms = 
tJ0.5(I xx + I yy ) and ellipticity ei). We measure the re- 
lease time r in units of "pixels," the amount of time it 
takes to clock the charge by a certain number of pixels 
in the parallel readout direction (i.e., one "pixel" is the 
time between row shifts during readout). 

In terms of astrometry, photometry, and size, there are 
two limiting regimes. Charge traps with very short re- 
lease times (or slow CCD readout) push charge from an 
object's leading edge onto its core, and drag core charge 
into a short tail. Both effects shift the object away from 
the readout register. The net effect also increases the ob- 
ject's size, because the core contains more charge than 
the wings. A small amount of flux can be lost from the 
wings into a tail, so Aflux is always slightly negative. 
However, the smoothing inherent in trailing correlates 
adjacent pixels and has the perverse effect of increas- 
ing the S/N. Note that the limiting behavior at low r 
is as expected: in our model, all traps inside a charge 
cloud carry an electron to the adjacent pixel at every 
clock cycle. In a real CCD, some charge may be released 



from very fast charge traps part- way through the 3-stage 
clocking cycle and returned to their original pixel. This 
process would lower the effective density of charge traps 
with low t. 

Charge traps with long release times (or fast CCD 
readout) steal flux primarily from an object's leading 
edge, and return it to the image in pixels well sepa- 
rated from the object. This stolen flux lowers the de- 
tection S/N. It also shifts the centroid as before, and 
decreases the size. For intermediate r, these effects are 
dominated by the addition of a tail, which increases the 
overall size. One curious dependency upon measurement 
method is that, while the rms size d rms decreases with 
r, the FWHM fitted by SExtractor increases: for exam- 
ple, AFWHM is negative for small r. This is presumably 
related to the net increase in detection S/N, and the seg- 
mentation of the image into fewer pixels that SExtractor 
determines belong to a given object. Note that if high r 
were achieved by dramatically speeding the CCD read- 
out, our assumption of instantaneous capture times may 
become invalid. A probabilistic capture mechanism over 
a finite time would result in lower effective densities of all 
traps, and potentially increased sensitivity to the density 
of charge throughout a pixel potential, changing the well 
filling parameters a and d. 

The spurious ellipticity induced in an object is interest- 
ingly different. The tail and the centroid shift induced 
by charge traps with short release times both elongate 
an object in the readout direction. As the tail lengthens, 
the spurious ellipticity initially increases. However, once 
the charge in the tail is sufficiently disconnected from 
the object and the object's centroid shifts back towards 
the correct position, the spurious ellipticity begins to de- 
crease. In the limiting case of charge traps with very long 
release times, charge missing from the object's leading 
edge could potentially elongate the object perpendicular 
to the readout direction; however, the residual centroid 
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Fig. 5. — The effect of charge trap release time r on a mea- 
surement of photometry, astrometry, size and ellipticity of a typ- 
ically small, faint galaxy degraded by CTI. Each y axis repre- 
sents the fractional change in that quantity. The absolute values 
of the y axes are largely irrelevant, depending upon the assumed 
density of charge traps, CCD well filling model, galaxy SN (13), 
size (FWHM=3.8 pixels) and morphology (circularly symmetric De 
Vaucouleurs profile). However, the trends reveal several lessons for 
future hardware: notice particularly the local maximum in |Aei|, 
which implies a worst case clocking time, or readout cadence, for 
CCDs. 

shift in this case is sufficient to maintain a small ellip- 
ticity in the readout direction (this result may depend 
upon the object's radial profile). 

Thus, we find that, in terms of weak lensing shear mea- 
surement, not all CTI is equally bad. Furthermore, there 
is a worst possible case, in which traps with release times 
corresponding to 3-4 clock cycles induce the most spuri- 
ous ellipticity. This value depends upon the shape mea- 
surement method: with KSB (Kaiser, Squires & Broad- 
hurst 1995) and RRG, it depends upon the size of the 
Gaussian weight function. The bump in Ad(r) around 
this value is real and also depends upon this scale. How- 
ever, from a more general argument about the dissocia- 
tion of flux from an object in a very extended trail, it is 
clear that a local maximum in |Aei| will be inevitable for 
all shear measurement methods. The clock speed is a pa- 
rameter that can be tuned in the hardware. We discuss 
this possibility in Sj6] 

5.2. Effects on galaxy morphology 

We create a series of images at the fiducial irradiation 
level of D08 (110 years at L2). We created galaxies with 
De Vaucouleurs (DvC) profiles and exponential profiles. 
For the DvC galaxies, we varied S/N, size, and input el- 
lipticity. For each different set of simulation parameters, 
we create 1000 simulated galaxies, each with a different, 
random sub-pixel position of the galaxy centroid and dif- 
ferent background noise realization. We measure the size 
in terms of the rms size d. Small galaxies have a size close 
to that of the PSF, representing the typical galaxies that 
will dominate a lensing survey; large galaxies are signif- 
icantly bigger than the PSF. The values of S/N, size, 
and ellipticity for the different simulations are shown in 
Table [5] The key result for weak lensing, the change in 
measured ellipticity, is illustrated in Figure [6] 

As expected, the small galaxies are significantly more 



affected by CTI than large galaxies. We also show that 
for small galaxies, brightness (S/N) is a mitigating fac- 
tor (but not for large galaxies); small bright objects are 
slightly less affected by CTI than small faint ones. An- 
other interesting feature recovered from these simulations 
is the dependence of |Ae| = |Aei| on galaxy ellipticity. 
Galaxies that are already aligned along the readout di- 
rection (ei < in the case of our simulations) are less 
affected by CTI than galaxies that are aligned perpendic- 
ular to the readout direction (ei > 0). Perturbations in 
the y direction (such as CTI) affect galaxies that are al- 
ready aligned in the y direction less than galaxies aligned 
in the x direction. 

5.3. Effects of trap density 

Figure also shows the results of simulations with 
varying trap densities. We increment the trap density (in 
all three species) from zero to the density that would be 
found after 220 years at L2 (twice the fiducial value from 
D08). At each of 50 evenly spaced points along this time- 
line, we create 100 simulated galaxies, each with a DvC 
profile (small, low S/N, e = 0). As expected, the degra- 
dation of shape measurement increases linearly with trap 
density. We find that 

= (2.65 ± 0.02) x 10" 4 [yr _1 at L2] (2) 

at 

for a radiation dosage averaged over an 11-year solar cy- 
cle and yearly displacement damage dose of 1.610 6 MeV 
per gram of Si. We make the simplifying assumption of 
zero traps (CTI= 0) at time t = 0. However, D08 show 
that the number of traps in a new LBNL CCD is smaller 
than the measurement error, so the approximation we 
make here is a good one. 

6. CONSEQUENCES FOR THE FUTURE SPACE MISSIONS 

Amara & Refregier (2008) ascertain that the multi- 
plicative error on measured shear needs to be kept below 
one part in 10^ 3 in order that future dark energy mis- 
sions not be dominated by systematic errors. This means 
that all sources of shear measurements error (not just 
the portion due to CTI) must be kept below Ae < 10~ 3 
throughout the mission lifetime. This level of shape mea- 
surement accuracy is represented in Figure |6] by a hori- 
zontal dotted line. 

The prediction for Ae after a fiducial 5 year mission 
can be calculated from Equation ©. To avoid any as- 
sumption of linearity with p, however, we have also run 
a larger number of simulations at 2 and 5 years of mean 
L2 exposure. In each case, we created 3000 simula- 
tions to reduce measurement noise due to the sub pixel 
galaxy position and sky noise. We find that Ae2 years = 
0.490 ± 0.01 x 10" 3 and Ae 5 years = 1.56 ± 0.02 x 10~ 3 . 
That is, without any correction, a 5 year weak lensing 
mission will have the entire shape measurement error 
budget consumed by CTI-induced effects before the end 
of the mission, even with specially designed, fully de- 
pleted, radiation hardened, p-channel CCDs. 

Fortunately, recent work using data from the HST's 
Advanced Camera for Surveys Wide Field Camera 
(ACS/WFC) has shown that software postprocessing 
can correct the effects of CTI on galaxy shapes by about 
a factor of 10 (Massey et al. 2009). Using the same 
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TABLE 2 

Parameters of Simulated Galaxies 



Profile 


S/N 


AIS/N) 


e 


Ae 


rms size d rrns [pixels] 


|A»| 


DvC 


low (13) 


l.l 





-0.028 


(small) 1.8 


0.19 


DvC 


low (13) 


l.i 


-0.17 


-0.024 


(small) 1.8 


0.19 


DvC 


low (13) 


1.2 


+0.17 


-0.031 


(small) 1.8 


0.19 


DvC 


high (50) 


0.7 





-0.024 


(small) 1.9 


0.17 


DvC 


low (20) 


0.3 





-0.010 


(large) 3.1 


0.20 


DvC 


high (50) 


3.0 





-0.010 


(large) 3.6 


0.20 


Exponential 


high (49) 


0.9 





-0.022 


(small) 1.9 


0.18 



0.008 



0.06 
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Fig. 6. — Left Panel: The degradation in shape measurements (|Aei|) as a function of time at L2. We make simulations for radiation 
exposure from to 220 years at L2 (the latter number being twice the accumulated radiation exposure used in D08). The data points 
with error bars are small, faint De Vaucoulcurs profiles far from the readout register, and therefore represent a worst-case scenario for 
CTE-induced shape errors. The other data points represent different levels of SN, size, and profile type as described in Table[2] We create 
and read out 3000 galaxies of each type, so that the error bars on those points are smaller than the plotting symbols. All of the individual 
profile types are plotted for 110 years of L2 radiation exposure; some points are slightly offset for readability. Right Panel: An enlarged 
portion of the left panel showing the effects of radiation damage in the first 20 years at L2. The data points at 2 and 5 years are described 
in J6]and have error bars about the size of the plotting symbols. The solid line is our best fit to our data points as given in Equation 2. The 
dashed horizontal line is at Ae = 0.001, the level of shape measurement accuracy needed for a future all-sky dark energy mission (Amara 
& Refregier 2008). 



code described here, trailed charge could be moved 
back to where it belonged in an iterative procedure 
restoring images to their true appearance. Massey et 
al. (2009) measured the time constants and number 
density of traps in ACS/WFC as a function of time, 
using extragalactic survey imaging that would be natu- 
rally available in any future survey without additional 
overhead. For this software, the factor of ten level of 
correction will be maintained down to the regime of 
future missions with much lower trap densities. One 
component of noise (due to variations in the number 
of traps in a given pixel, which we treat as a constant 
density) will be improved, but this affects only scatter 
in Ae rather than the level itself. We therefore conclude 
that, using CCDs with the characteristics of those in 
our study, and proven software mitigation techniques 
to achieve an additional factor of 10 correction, CTI 
in a future dark energy mission would be satisfacto- 
rily controlled at only 10% of the total shape error budget. 

One way in which CTI models (and mitigation tech- 
niques) could potentially be pushed beyond correction 
by a factor of 10 would be to precisely locate individ- 



ual charge traps, rather than treating them statistically. 
This would be most beneficial in the early years of a dark 
energy mission, when p « 1. Designing flexibility into 
the clocking speed, waveform, and voltage in CCD elec- 
tronics provides the ability to locate traps via "pocket 
pumping" (Janesick 2001). Pocket pumping is a process 
in which a uniform level of charge introduced by a flat 
field lamp is rapidly shuffled back and forth thousands 
of times in the parallel direction. Following the charge 
shuffle, the charge is transferred to the readout transistor 
in the normal manner. The resulting image reveals accu- 
mulated charge captured and then released by each trap 
in the shape of a dipole of an overdensity neighboring an 
underdensity of charge. The orientation and strength of 
this dipole reveal the location of the trap within the pixel 
and the effectiveness of the trap. This can be repeated 
with different levels of initial charge to map out traps 
in 3 dimensions within the CCD. Because the readout 
time of the CCD is dominated by the clocking of serial 
charge, the pumping of charge in the parallel direction 
does not introduce a significant amount of overhead to 
the survey. For example, it takes approximately 45 sec- 
onds to record a normal image from the LBNL CCDs 
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using the D08 clocking parameters. Image acquisition 
takes an additional 45 seconds for 20000 cycles of pocket 
pumping with a five pixel shift. Acquisition of five succes- 
sive pocket-pumping images every week to reject cosmic 
rays and model the traps therefore would account for less 
than ten minutes of additional calibration time. 

Another suggested mitigation technique is "charge in- 
jection" (also called "fat zero" or "pre-flash") in which 
charge is placed into the pixels in order to fill the traps. 
The charge injection is simply a method to increase the 
overall background level and fill the volume. However, 
this has the effect of increasing the background in an im- 
age (much like increasing the zodiacal background) and 
will reduce the S/N of the detected objects. This is obvi- 
ously undesirable for a weak lensing experiment in which 
the observer is attempting to measure shapes of faint 
galaxies. 

We showed in §5. II that there exist clocking time scales 
that are maximally bad for shape measurement. If the 
clock cycle is 3-4 times the trap release time, then Ae is 
maximized. Thus, future missions with weak lensing as 
a primary science driver should include an optimization 
of the charge clocking time in their CCD readout elec- 
tronics. Increasing the rate at which charge is clocked 
serially, and thus increasing the rate at which parallel 
transfers of rows can be made, increases readout noise, 
resulting in an effective loss of survey depth; decreas- 
ing the rate at which charge is serially clocked increases 
the readout time and, if that dominates over factors like 
slew and settle time between exposures, will reduce sur- 
vey area. Thus, careful consideration must be paid to 
the trade-offs in any such optimization. Furthermore, 
we are only making recommendations for how to mini- 
mize shape measurement errors due to CTI. Photometry 
is also degraded by traps with large values of r relative 
to the charge clock period. This can remove charge from 
objects but place it far enough away that the shape is not 
significantly affected. Thus, trade-offs in charge clocking 
time must also take into account the photometric accu- 
racy requirements of the mission. 

We finally note that the temperature at which detec- 
tors are operated at has significant effect on CTE, and 
thus future missions should be tested and optimized with 
this in mind. 

7. CONCLUSIONS 

We have quantified the effect of CTI on measurements 
of weak gravitational lensing. We first simulated the 
transfer of charge within LBNL p-channel CCDs using 
a model for charge traps with three characteristic re- 
lease times to reproduce the experimental results of D08. 
Using this model, we then simulated deep exposures of 
galaxies for weak lensing measurements from a space- 
based telescope subject to radiation damage. The re- 
sulting simulated data were used to quantify the effects 
of radiation damage on shape measurements of galaxies 
of various sizes and S/N levels, the true shapes of which 
were precisely known. Most galaxies in any weak lensing 
survey will be small and faint; as expected, we find that 
these suffer worst from the effects of CTI. 

The level of CTI-induced shape error Ae will approach 
the total shape error budget of a dark energy mission (1 
part in 10 -3 ; Amara & Refregier 2008) after less than 4 
years of radiation exposure at L2. Software mitigation 



techniques in image postprocessing, already proven on 
HST data (Massey et al. 2009), will be able to reduce the 
levels of shape error well below mission requirements. 
We have also suggested hardware capabilities, such 
as "pocket pumping" and adjustments to the readout 
speed, that may provide additional help. However, our 
numerical results are only valid for p-channel devices, 
whose CTE characteristics after radiation exposure have 
been shown to be superior to more common n-channcl 
devices (Lumb 2009; Marshall et al. 2004). Given the 
necessity of both hardware and software mitigation of 
CTI effects for successful mission operation, we recom- 
mend that future spacecraft be designed with detectors 
and mission parameters that ensure CTE characteristics 
no worse than the LBNL p-channel devices we simulated 
for this work. 

There are several caveats to our results. First, we as- 
sumed that all galaxies are small, faint, and lie far from 
the readout register. Galaxies that are bright, large, or 
nearer the readout register will suffer less from CTI. In- 
deed, the average distance to the readout register will 
be exactly half of the worst case scenario we outlined, 
so the mean Ae will be a factor of two lower. In a real 
mission, there will be several dithered exposures of each 
galaxy, with each dither placing the galaxy a different 
distance from the readout register; this may allow us to 
further model the effects of CTI on shapes and partially 
mitigate those effects. We also assumed a high level of 
radiation from the Sun for at least a portion of the mis- 
sion due to modeling radiation flux in the heaviest part 
of a typical 11-year solar cycle. A mission flown during 
the level of minimum particle radiation at L2 would suf- 
fer less radiation damage, but a mission flown entirely 
during the maximum of the solar radiation flux would 
suffer more damage. This is a large uncertainty because 
the radiation flux due to the Sun can vary by an order 
of magnitude over the solar cycle (Barth et al. 2000). 
Furthermore, any real mission would need to adjust the 
flux according to the planned shielding on the space- 
craft (D08 assume the SNAP design) and should take 
into account secondary particle cascades from reactions 
of high energy radiation with the shielding material (D08 
ignored such secondary radiation). The 'stacking' and 
'trails' methods probe different properties of CTE and 
we have chosen the fits to trails in D08 to model the ef- 
fects of radiation exposure on shape measurements. We 
have only explored the charge transfer and radiation tol- 
erance properties of a certain model of CCD operating at 
a single temperature; any future space-based weak lens- 
ing missions should undertake a similar analysis using 
the CCDs planned for that mission. Finally, if the CCDs 
contain a significant density of traps even before they are 
launched into the harsh radiation environment of space, 
the CTI will be worse throughout, and the useful mission 
lifetime reduced. Clearly, then, this paper is a first step, 
and any future mission should use a procedure similar 
to the one we have developed in this paper to meet spe- 
cific mission requirements by fully optimizing its choice 
of CCDs, clocking rate, and shielding. 
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